LOCAL CONSERVATION LAWS AND THE HAMILTONIAN 
FORMALISM FOR THE TODA HIERARCHY REVISITED 



FRITZ GESZTESY AND HELGE HOLDEN 



Abstract. We revisit an elementary recursive approach to local conservation 
laws and the Hamiltonian formalism for the Toda hierarchy. 



1. Introduction 

In this paper we revisit local conservation laws and the Hamiltonian formalism for 
the Toda hierarchy. Consider sequences {a(n, to), b(n, to)} n ez £ £°°(Z) satisfying 
some decay assumptions as n — » ±00 to be specified later, parametrized by the 
deformation (time) parameter to £ R, that are solutions of the Toda equations 



v 6 t0 -2(a 2 -(a-) 2 ) 
Here denote shifts, that is, c ± (n) = c(n ± 1), n G Z. Then clearly 

d t0 ln(o(*o, n)) = 9 t0 £ 6(t , n) = 0, (1.2) 



provided one has sufficient decay of the sequences a, b. Additional calculations also 
show that 



dt ( a (*o, ™) 2 + \b(t , n)A = 0. (1.3) 



Indeed, one can show the existence of a sequence {p+j jjgN of polynomials of 
a, b and certain shifts thereof, with the property that the lattice sum is time- 
independent, 

d to J2 P+>j (t ,n)=0, jeN. (1.4) 

nez 

This result is obtained by deriving local conservation laws of the type 

d toP+tj + (S + -I)J + , 0d = 0, jeN, (1.5) 

for certain polynomials J+,o,j of a, b and certain shifts thereof. The polynomials 
J+,o,j constructed via a straightforward and explicit recursion relation. 
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The above analysis extends to the full Toda hierarchy as follows. The pth equa- 
tion, p 6 No = N U {0}, in the Toda hierarchy is given by 

T\Aa,b)=h~ aifp++1 ~ U+l) )=0, t p eR, (1.6) 

\ K + g P +i - g p+ i J 

where fi, ge are carefully designed polynomial expressions of a, b and certain shifts 
thereof. Recursively, they are given by 

fo = 1, ,g = -ci, (1.7) 
2f i+ i +9i + 9J -26/* = 0, feN„, (1.8) 
g i+ i-9i +1 + 2(a 2 f + -(a-) 2 fl)-b( gi -g^) = 0, £eN Q , (1.9) 

where c± S C is a constant. On each level in the recursion an arbitrary constant 
C( G C is introduced. In the homogenous case, where all these constants ci, £ € N, 
are set equal to zero, a hat is added in the notation, that is, fi,gi, etc., denote 
the corresponding homogeneous quantities. The homogeneous coefficients fe,ge 
can also be expressed explicitly in terms of appropriate matrix elements of powers 
of the Lax difference expression L defined in (|1 . 1 If) below (cf . Lemma I2.2|l . The 
conserved densities p+.j are independent of the equation in the hierarchy while the 
currents J+, P ,j depend on p; thus one finds (cf. Theorem 14. 6|) 

9t P P+,j + (S+ - I)J+, PJ =0, j e N, p e N . (1.10) 

Our approach is based on a careful analysis of the high-energy expansion of the 
Green's function for the operator H in i? 2 (Z) corresponding to the difference ex- 
pression 

L = aS+ + a-S- +b. (1.11) 

In particular, we do not depend on real-valued sequences a, 6, and thus the operator 
H is permitted to be non-self-adjoint. 

In addition, we study the Hamiltonian formalism for the Toda hierachy. Proving 
the relations (cf. Lemma 

^ = --&-!, (1.12) 
da a 

y±=£f t _ lt leN, (1.13) 
one infers that the pth equation in the Tode hierarchy can be written as (cf. Theorem 

(here N_i = NU{0,-1}), where the functional I p are given by 

p 1 
I-i=I-i, I p = ^2 c p -kh, P £ No, I q = —^^2 (f q+2 {n) - \ q+2 ), q eN-i. 

k=0 1 neZ 

(1.15) 

Here D is an explicit 2x2 matrix- valued difference expression (cf. (|5.7ll ) and A p+ 2 
are counterterms necessary to guarantee convergence of the infinite series involved. 
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Furthermore, we show that any I r is conserved by the Hamiltonian flows in Ijl.l4|l (cf. 
Theorem 15 .6|) . that is, 



Moreover, for general sequences a, b ( i.e., not assuming that they satisfy an equation 
in the Toda hierarchy), we show in Theorem 15. 71 that 



for a suitably defined Poisson bracket { • , • } (see ()5.f 9(1 ). 

Finally, we also show the existence of a second Hamiltonian structure for the 
Toda hierarchy (cf. I(5.55|l . l(5.56jl L and sketch the possible extension of this formal- 
ism to the case of almost periodic sequences a, b on Z. 

The Toda hierarchy and its interpretation as a Hamiltonian system has of course 
received enormous attention in the past. Without being able to provide a com- 
prehensive review of the vast literature we mention that Lie algebraic approaches 
to Hamiltonian and gradient structures in the asymmetric nonperiodic Toda flows, 
Poisson maps, etc., are studied, for instance, in 0]. Master symmetries for the 
finite nonperiodic TL are treated in [8], the 2nd Hamiltonian structure for the 
periodic Toda lattice is discussed in [3], and multiple Hamiltonian structures for 
various Bogoyavlensky-Toda systems, Poisson manifolds, master symmetries, re- 
cursion operators, etc., can be found in [Hj and ^jj. Moreover, Poisson brackets are 
studied in ^3 Sect. III. 4], |17| . |27l Sect. II. 4], and constants of motion appeared 
for instance, in O Sect. 7.3], Q3J Sect. III.2.4], US], [2U|, EH Sect. 13.4], and [321 
Sect. 3.7]. Fairly recently, an elementary approach to the infinite sequence of Toda 
hierarchy conservation laws was presented in |37j (but it seems to lack an explicit 
and recursive structure). Moreover, a Hamiltonian formalism for the Toda hier- 
archy, based on recursion operators and hence on an algebraic formalism familiar 
from the one involving formal pseudodifferential expressions in connection with the 
Gelfand-Dickey hierarchy, was recently developed in [35], 

Still, we emphasize that our recursive and most elementary approach to local 
conservation laws of the infinite Toda hierarchy does not appear to exist in the 
literature. Moreover, our treatment of Poisson brackets and variational derivatives, 
and their connections with the diagonal Green's function of the underlying Lax 
operator, now puts the Toda hierarchy on precisely the same level as the KdV 
hierarchy with respect to this particular aspect of the Hamiltonian formalism. 



In this section we briefly review the recursive construction of the Toda hierarchy 
and associated Lax pairs following 0, [221 Sect. 1.2], and O Ch. 12]. 
Throughout this section we make the following assumption: 

Hypothesis 2.1. Suppose 




(1.16) 



{I p J r } = 0, p,r£ N_i, 



(1.17) 



2. The Toda hierarchy in a nutshell 



a,b£l°°(Z), a(n) ^ for all n £ Z. 
We consider the second-order Jacobi difference expression 



(2.1) 



L = aS + +a~S~ +6, 
where S denote the shift operators 

(S ± f)(n) = /="=(») = f(n±l), n £ Z, / £ C z . 



(2.2) 



(2.3) 
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Here C J denotes the set of complex- valued sequences indexed by J C Z. 

To construct the stationary Toda hierarchy we need a second difference expres- 
sion of order 2p + 2, p G No, defined recursively in the following. We take the 
quickest route to the construction of P 2p+2 , and hence to the Toda hierarchy, by 
starting from the recursion relations (|2.4|) - l|2.6[) below. 

Define {fe}een and {gij^tio recursively by 

fa = 1, 9o = -ci, (2.4) 
2ft+i+9t+97-2bf t = 0, feNo, (2.5) 
ge+i-gi +1 + 2(a 2 f+ -(a-) 2 f e -) -b(g e ~gi) = 0, I e N . (2.6) 
Explicitly, one finds 

/o = l, 
fi =b + cx, 

f 2 = a 2 + (a-) 2 + b 2 + Cl b + c 2l etc., (2.7) 

go = -ci, 

gi = -2a 2 - c 2 , 

g 2 = ~2a 2 (b + + b) + ci(-2a 2 ) - c 3 , etc. 

Here {cj}jeN denote undetermined summation constants which naturally arise when 
solving $Tty - $nty . 

Subsequently, it will be convenient to also introduce the corresponding homoge- 
neous coefficients fi and gi, defined by vanishing of the constants fc£N, 

/° = 1 > ft = Mc>=a,k=i,...,t> ^ eN ' 

& = «L=o,k=w+i> £eN °- ( 2 - 8 ) 



Hence, 



introducing 



ft = ^2 c ?-kfk, 9i = ^ Q-fcfffc - Q+i, I £ N , (2.9) 



fc=0 fc=i 



c = 1. (2.10) 
The following explicit characterization of fi, gt turns out to be useful. 

Lemma 2.2. Suppose a,b satisfy Hvvothesis \2.1l Then the homogeneous coeffi- 
cients {fi}ee® and {gt}l&i satisfy 

fi{n) = (5 n ,L e 8 n ) eHZ) , £eN , (2.11) 

g e (n) = -2a(n)(S n+1 ,L e S n ) e 2 iz) , I e N , (2.12) 

where S n = {<5„ jm } me z, n6 2, and ( • , • )^2( Z ) denotes the usual scalar product in 

e 2 (z). 

For the proof see [5] . 

Next we define difference expressions P 2p + 2 of order 2p + 2 by 

p 

P 2p+2 = -LP +1 + Y / (gi + '2af e S + y p - e + f p+1 , P eN . (2.13) 
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Introducing the corresponding homogeneous difference expressions P 2p +2 defined 

by 

P21+2 = P2i+2\ Ck=o k=1 t e No, (2.14) 

one finds 

p 

P2 P +2 = c p -£P 2 e+2- (2-15) 
1=0 

Using the recursion relations 12.4Jl - H2.6|l . the commutator of P 2p + 2 and L can be 
explicitly computed and one obtains 

[P 2p+2 , L] = - a{ 9 ; + g p + f+ +1 + f p+1 - 2b+f+)S+ 

+ 2( - b(g p + f p+1 ) + a 2 f+ - {a-ff- + b 2 f p ) 

-a-(g p + g- +f p+1 + f- +1 -2bf p )S-, p 6 No. (2.16) 

In particular, (L, P 2p + 2 ) represents the celebrated Lax pair of the Toda hierarchy. 

Remark 2.3. Varying p 6 No, the stationary Toda hierarchy is then defined in 
terms of the vanishing of the commutator of P 2p +2 and L in i|2.16[l by 

[P 2p+2 ,L] = s-TL>, b) = 0, pe N . (2.17) 

Equation H2.17[l can be shown to be equivalent to 

s-Tl p (a, b) = ( f P +1 ~ fp _ +1 ) = 0, pe N . (2.18) 



9 P +i ~ 9 P +i 



Specifically, one obtains 



s-Tlo(a, 6 )=( 2((a 6 : )2 _ 5 a2) )=0, 

sTUab)- ( (* + ) 2 -(--) 2 + (b + ) 2 -b 2 \ (21Q) 
s lb(a, b) - ^ 2(a -)2( 6 + 6 _j _ 2a 2 {b + +b) j (2.1Jj 

+ Cl (2((a-) 2 "-« 2 )) =0, 6tC -' 
for the first few equations of the stationary Toda hierarchy. 

Next, we introduce polynomials F p (z,n) and G p +\{z,n) of degree p and p + 1 
with respect to the spectral parameter z G C by 

p 

F p (z,n)=Y.fp-t( n > e > ( 2 - 2 °) 
£=0 

p 

G p+ i{z, n) = - Z p +1 + J2 9p-eHz e + f P+1 (n). (2.21) 

Explicitly, one obtains 
Fo = 1, 

Fx = z + b + c x , 

F 2 =z 2 + bz + a 2 + (a~) 2 +b 2 + Cl (z + b)+c 2 , etc., (2.22) 
G x = -z + b, 

G 2 = -z 2 + (aT) 2 - a 2 + b 2 + cx(-z + 6), etc. 
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Next, we study the restriction of the difference expression P2 P +2 to the two- 
dimensional kernel (i.e., the formal null space in an algebraic sense as opposed to 
the functional analytic one) of (L — z). More precisely, let 

ker(L -z) = {y;:Z^CU {00} \ (L - z)iP = 0}. (2.23) 

Then (f2~T3|) implies 

P2 P+ 2 |kcr(L- 2 )= {2aF p (z)S+ + G p+1 (z))\ kci{L _ z) . (2.24) 

Next we turn to the time-dependent Toda hierarchy. For that purpose the func- 
tions a and b are now considered as functions of both the lattice point and time. 
For each equation in the hierarchy, that is, for each p, we introduce a deforma- 
tion (time) parameter t p £ R in a, 6, replacing <z( • ), b( • ) by a( ■ , t p ), b( ■ , t p ). The 
second-order difference expression L (cf. (12.2(1 ) now reads 

L(tp) =a{-,t p )S + + a~(-, tp)S~ +b(-,t p ). (2.25) 

The quantities {fe}i£N , {ge}een , and P2 P +2, P £ No are still defined by (|2.4(l - l|2.6(l 
and l(2.13|l . respectively. The time-dependent Toda hierarchy is then obtained by 
imposing the Lax commutator equations 

L tp (t p )-[P 2p+2 {t P ),L(t p )} = 0, t p e R, (2.26) 

varying p £ No- Relation 12.26|l implies 

(a tp + a(g+ + g p + f+ +1 + f p+1 - 26+/+)) S+ 

- (-b tp + 2( - b(g p + f p+1 ) + a 2 f+ - (a-ff- + b 2 f p )) (2.27) 

+ (a tp + a( 9 ; + g p + f+ +1 + f p+1 - 2b+f+))~ S~ = 0. 

In addition one computes 

= L tp -[P 2p+2 ,L] 

= (a tp - a(/+ +1 - f p+1 )) S + - (-b tp - g p+ i + g^) 

+ K ~ a(f+ +1 - f p+1 ))- S- . (2.28) 

Varying p G No, the collection of evolution equations 

Tl p (a, b) = \ tp ~ a(/p+1 ~ /p+l) ) = 0, (n, t p ) e Z x R. p e N (2.29) 
V K + 9 P +i - 9 P+ i J 

then defines the time-dependent Toda hierarchy. Explicitly, 



Tl (a,&) - 
Tli (a, b) = 




(2.30) 



represent the first few equations of the time-dependent Toda hierarchy. The system 
of equations, Tlo(a, 6) = 0, is of course the Toda system. 

The corresponding homogeneous Toda equations obtained by taking all summa- 
tion constants equal to zero, ct = 0, £ = 1, . . . , p, are then denoted by 

T\ p (a,b)=Tl p (a,b)\ c ^ oi=l tp . (2.31) 
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Restricting the Lax relation (|2.26(l to the kernel ker(L — z) one finds that 

' lker(i- Z ) = ( Lt p + ( L - Z ) P 2p+ 2 ) \ kcr{L _ z) 

2(z - b+)F+ - 2{z - b)F p + G+ +1 - G- + ^j S + 
+ (b tp + (z- b)^ + 2(a-fF p - 2a 2 F+ + (z - b)(G p+l - G p+l ) 



0=(L tp -[P 2p+2 ,L 



Hence one obtains 



a 



A = - 2 {z - b+)F+ + 2{z - b)F p + G~ +l - G+ +1 , 

<z _ & )% _ 2(a-fF- + 2a 2 F+ - (z - b)(G p+1 - G p+1 ) 



Further manipulations then yield, 

a tp = -a(2(z - b+)F+ 



G p+ i + G p+ i) 



b tp = 2((z - bfF p + (z- b)G p+1 + a 2 F+ - (0 2 F"). 



(2.32) 



kcr(L — z) 

(2.33) 
(2.34) 

(2.35) 
(2.36) 

Indeed, (|2.35(1 follows by adding G p +i — G p +i to (|2.33(1 (neglecting a trivial summa- 
tion constant), and an insertion of (|2.35(l into (|2.34(1 implies 1(2.361) . Varying p e No, 
equations (|2.35|) and ((2.361) provide an alternative description of the time-dependent 
Toda hierarchy. 

Remark 2.4. From ((2.4(1 - ((2. 6|) and (J220J), (12.21(1 one concludes that the coeffi- 
cient a enters quadratically in F p and G p +\, and hence the Toda hierarchy ((2.29(1 
[respectively 1(2. 18(1 ) is invariant under the substitution 

a £ = {e{n)a{n)} n& , e(n) S {1, -1}, n 6 Z. (2.37) 

3. Green's functions and high-energy expansions 

In the following it will be necessary to associate certain operators to the dif- 
ference expression L in 12.2(1 . Throughout this section we assume that a, b satisfy 
Hypothesis 12.11 

We start with the Jacobi operator H defined on £ 2 (Z) by 

Hf = Lf, /edom(if)=l 2 (Z). (3.1) 

In addition, we introduce the half-line Dirichlet operators on £ 2 ([no, oo) flZ) 

and no on £ 2 ([n Q , -oo) l~l Z), n e Z, by 



_ J a(n )u(n + 1) + b(n )u(n ), n = n 0} 

1 (Lu)(n), n E [no + 1, oo) fl Z, 

u e dom(Hf nQ ) = £ 2 ([n , oo) n Z), (3.2) 

a(n )v(n a + 1) + b(n a )v(n a ), n = n , 

(Lv)(n), n G (— oo, no — 1] fl Z, 

v e dom(H D no ) = £ 2 ((-oo,n a ] nZ). (3.3) 

(It is customary to join the Dirichlet-type boundary condition u(no — 1) = in the 



(H° no u)(n) 



(HB ino v)(n) = 



case of H? n , and v{uq + 1) = in the case of H_ n in order to formally be able 
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to avoid the case distinction in l|3.2|) and (|3.3Jl . s ) The following elementary result 
about a spectral inclusion of H and H± n will be useful later. 

Theorem 3.1. Suppose a,b satisfy Hvvothesis \2.1\ Then the numerical range of 
H , and hence in particular, the spectrum of H and H± n is contained in the closed 
ball centered at the origin of radius 2 1| ci, || oo + ||6||oo, that is, 



spec(H), spec(fr£ no ) C B(0; 21101100 + (H^), n G Z. (3.4) 

Proof. We denote by W(T) the numerical range of a bounded linear operator T G 
B(H) in the complex, separable Hubert space H, 

W(T) = {(g, Tg) H G C | 5 G \\g\\ n = 1}. (3.5) 

It is well-known that W(T) is convex and that its closure contains the spectrum of 
T, 



spec(T) C W(T) (3.6) 
(see, e.g., [231 Ch. 1]). Elementary arguments then prove 

l(/,^/)^(z)l<(2||o||oo + ||6||oo)||/|||i (Z ), fe£ 2 (Z), (3.7) 

\(w,H^ no w) £2([no±oo)nZ) \ < (2||o||oo + ||6||oo)||w||?.( [no ,± 00 )n Z ). (3-8) 

w G ^ 2 ([n ,±oo) HZ), 

and hence IpTl)) . □ 

Since if is a bounded second-order difference operator in £ 2 (Z) with 

||-ff|| < 2HIOO + llfcHoo, (3.9) 

the resolvent (H — zl)^ 1 , z G C\spec(_ff), of H is thus a Carleman integral operator 
(cf., e.g., Sect. 6.2]; the corresponding measures involved are of course discrete 
measures). In particular, introducing 

ip+(z,n) = {{H - ziy^a)^), n>l, z e C \ spec(iJ), (3.10) 

i>-(z,n) = {(H - ziy^o)^), n<-l, z e C \ spec(iJ), (3.11) 

where 

;iez\ W , fcez ' (3 - i2) 

and using the Jacobi equation Lip±(z) = zip(z) to uniquely extend ^p±{z,n) to all 
n G [0, Too) n Z (this is possible since by Hypothesis 12. II a(n) ^ for all n G Z), 
one obtains, 

Ltp±(z) = zip±(z), ip±(z, •) G ^ 2 ([n ,±oo) HZ), z G C \ spec(-ff), n G Z. 

(3-13) 

We note that the Weyl-Titchmarsh-type solutions ip±(z, ■ ) are unique up to nor- 
malization. Moreover, by the second inclusion in l|3.4|) . 

i>±(z, ■ ) is zero free for \z\ > {2\\a\\oo + \\b\\oo) (3.14) 
since ip±(z, no) = for some no G Z yields an eigenvalue of H± no ±i- 
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The Green's function G(z, • , • ) of H, that is, the integral kernel of the resolvent 
(H — zT)^ 1 (with respect to discrete measures), is then given by 



G(z, to, n) 



G_„ (z,m,n) 



1 J ip-(z, m)ip+(z,n), m < n, 

W(ip-(z),t(j + (z)) )i/)-(z, n)ip + (z,m), m>n, (3.15) 
z G C \ spec(ff), to, neZ, 

where W(f, g)(ri) denotes the Wronskian 

W(f,g)(n) = a(n)(f(n)g(n+l)-f(n+l)g(n)), n G Z, (3.16) 

of complex- valued sequences / = {f(n)} ne z and g — {g(n)} ne z C C. 

The corresponding Green's functions G± no (z, ■ , • ) of H± nQ , that is, the integral 
kernels of the resolvents (H — z/) _1 (again with respect to discrete measures), are 
then given by 

G D (zmn) = — [(po(z,m,n - l)ip+(z,n), m<n, 

+ '™° ' ' a(n Q -l)ip + (z,n -l))(f) (z,n,no-l)ilj + (z,m), m>n, 

z G C \ spec(if£ no ), to, n G Z, (3.17) 

— 1 J ip-(z, m)(j)o{z, n, no + 1), m < n, 

a(no + l)ip-(z,no + 1) 1 ^-(z, n)cf)o(z, to, n + 1), to > n, 

z G C \ spcc(H D ng ), to, 7i G Z, (3.18) 

where 4>o(z, ■ , no) is a solution of 

LV = *V>, (3-19) 

satisfying the initial conditions 

(f>o(z,n ,no) = 0, 0o(z, n + l,n ) = 1, z G C. (3.20) 
Next, assuming that -0 satisfies (|3.19f) . the function <p = <f>{z, n), introduced by 

4>=—: (3-21) 

satisfies the Riccati-type equation 

ac/j + a-^-)- 1 + (b- z) = 0. (3.22) 

Defining 

<f>±(z,n) = ^±( z ' n+1) , zeC\spcc(^„ ±1 ),neZ, (3.23) 
ip±{z,n) 

and using the Green's function representations l|3.17|l . (|3.18|) . one computes 
c/) + (z,n) = -a(n)G% !n+1 (z,n + l,n+ 1) 

= -a(n)(<W, (Jf£ n+1 - ^)- 1 )^+i), 2([n+1>oo)nz)! (3-24) 
2 G C \ spec(iJ^„ +1 ), n G Z, 

and 

</>_ (z, n) = a(n) _1 (z — 6(n) + a(n — 1) 2 G^ n+1 (z, n — 1, n — 1)) 

= a(nT 1 (z-&(n) (3.25) 
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+ a(n-l) 2 (6 n ^ 1 ,(H^ n _ 1 -zI)- 1 6 n ^) e2{{ _ oon _ 1]nz ^, 

z G C\spec(i/ D „_ 1 ), neZ. 

We provide the following description of the asymptotic behavior of <p± as z — > oo. 

Lemma 3.2. Suppose a,b satisfy Huvothesis \2.1\ Then 4>± has the following con- 
vergent expansion with respect to 1/z around 1/z = 0, 



7=1 



'+,3 



-3 



<M*) = < ; _, (3-26) 



where 



<t>+,i = 1, 0+,2=& + , 

3-1 

= &+</>+ j + (a + ) 2 E <^-^+,*> J ^ 2 ' ( 3 ' 2? ) 

£=1 

0- ,-i = 1, 0-.o = -b, <f>-,i = -(a - ) 2 , 

3 

= ~ H <t>-j-rf-4> 3 ^ L ( 3 - 28 ) 

Proof. Since if± n±l are bounded operators, (|3.24|) and (|3.25|) prove the existence of 
an analytic expansion of </>±(z) with respect to 1/z around 1/z = for |z| sufficiently 
large. Moreover, H3.24[l and (|3.25(1 yield the leading asymptotic behavior, 

<j>+(z) = - az- 1 + 0(z~ 2 ) (3.29) 

z — >oc 

and 

0_(z) = - + 0(1). (3.30) 

z^oo a 

Thus, by making the ansatz (|3.2tj|) for <f>± and inserting it into l|3.22|l one finds 
(|X27jl and JS^S). □ 

For the record we note the following explicit expressions, 

</>+,i = !> 
0+,2 = 6 + , 

+ , 3 = (a+) 2 + (6+) 2 , 

<A+,4 = 

<t>-,-\ = 1, (3.31) 
0-,o = -b, 
cb-,i = -(a-) 2 , 



(a+) 2 (26+ + 6++) + (6+) 3 , etc. 



,2 = -(0 2 &~, 

,3 = -(a-) 2 ((a-) 2 + (6-) 2 ), etc. 
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Later on we will also need the convergent expansion of ln(0+) with respect to 
1/z for l/\z\ suficiently small and hence we note that 



lnO + (z)) = In ( a^2<j) +J z- 
V 3=1 

= In (l \ + In + 4>+, ]+ iz 



- z , 
where 

= 0+,2, = <A+,j+i - V, ~ <ft+,i+i-lP+,i, j > 2. (3.33) 

£=1 ^ 

The first few explicitly read 

P+,i = b + , 

P + ,2 = {a + ) 2 + \{b + )\ (3.34) 
P+,s = l(b+) 3 + (a+) 2 (b+ + b++), etc.. 



Similarly, one finds for l/\z\ sufficiently small, 

, oo 

M<f)_(z)) = In (cT 1 J2 0-, 



j=-i 

f 7\ 

In 



i)+ln(l + ^0_ li _ 1 *-i) 



ln (-J +E^ z ~ J < ( 3 - 35 ) 



where 



j-i - ~j 0-,j-i-^p-^, i ^ 2 - ( 3 - 36 ) 



The first few explicitly read 

P-,l = -b, 

P-,2 = -(a- f - \b 2 , (3.37) 
P - :3 = -^b 3 -{a-) 2 {b- +b), etc. 
Equations l|3.33|l and (|3.3t)|) may be written as 

P±,i = 0±,i±i, P±,j = 0±,i±i ~ ~ <P±,j±i-£P±A J > 2 - (3.38) 
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The next result shows that fj and ±jp±j are equal up to terms that are total 
differences, that is, are of the form (S f+ — I)d±j for some sequence d±,j. The exact 
form of d±.j will not be needed later. 

Lemma 3.3. Suvvose a.b satisfy Hvvothesis WH\ Then, 

fj = ±jp±j + (S+-I)d±j, j6N, (3.39) 

for some polynomials d±j, j S N, in a and b and certain shifts thereof. 

Proof. To shorten the notation we introduce the abbreviations 

g(z, n) = G{z, n, n), h(z, n) = G(z,n, n + 1), (z,n)GCxZ. (3.40) 

Thus, 

J>+ = -• (3.41) 
5 

To increase readability we suppress the display of variables in the subsequent com- 
putations. First, one observes that 

a 5ff -MaA-l)=a ^ W^T - ^ 

= (a^+^i - V- + - (3.42) 



and thus. 



This yields 



h 2 ah 2 g + ah g + 



t 2 n/m+ 



agg+ g ah - 1 g 



(3.43) 



1 d f ah \ Id, fg+\ 

— ln(6 + ) = ~ — In — +-— In — . 3.44) 

dz w+y 2dz \ah-l) 2 dz \ g ) y ' 

Next we study the first term on the right-hand side of 13.44fl . namely, 

d ( ah \ h z ah z 

Tz ln {^h—i) = T-^h—Y (3 - 45) 



' - (1 - 2ah) + 4 ) + 4ah t . (3.46) 



We claim that 



h ah — 1 \ g g 

To this end we first take the logarithmic derivative of 1)3.42(1 . that is, of agg 
h(ah — 1), to obtain 

9z , gf h z ah 



Thus, 



g g + h ah — 1 



(3.47) 



h ah — 1 \ \ g <?+ 



h ah — 1 \ ft aft, — 1 

-- ^-(l + (a/i-l)-a/i) =0. (3.48) 
ah — 1 v 7 
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Combining 44j) . (|3.45|1 . and (|3.46j) one finds 

Considering the first term on the right-hand side of i|3.49[) and subtracting g then 
yields 

I . . ■ ) (9 1 + 9_l 



~(l-2afc) 2i + 2£- )+2ah z -g 
2 \9 9 + J 

= ~(1 - 2a/i)f — ln(a/i + a~/r - 1) + — \n{a + h + + ah - 1) 
2 \ dz dz 

d , ,. ,.,\ _ , ah + a~h~ — 1 



ln(0-&+)(2-&))^ + 2a/i 



dz J z — b 

1 (1 2ah)( lr - " i ' 



2 \ah + a h —1 a+/i+ + a/i — 1 V z ~ ^ + z — 6 

ah + a~h~ 1 
2a/i, 



z — 6 z — & 
1 \ o/i a _ /i" 



2 \ z — & z — 6+ / z — 6+ z — & 

„ , 1/, „ , % / ah z + a~h~ a + hf+ah z \ ,„ „„. 

+ 2a/i* + - 1 -2a/i — — 2— + — — ^ , 3.50 

2 V ; V a/l + a ~^~ ~ 1 a+/i+ + a/i-l/ v ' 

using 

- % - aft - aT/T + 1 = 0. (3.51) 

The latter follows from the definition l|3.40|l of g and h (cf. Q3.15I) - ) and from Lt/>_ = 
zip—. By purely algebraic manipulations one obtains that 

1, v / ah T -\-arh~I a + hf + ah 

2ah z + -(l-2ah)' 



2 \ah + a-h--l a+h + + ah-l 

a~h~ah z — aThZah aha + hf — ah z a + h + 



ah + a h — 1 a + h + + ah — 1 

1/ a + ht — ah z ah z — a~h~ 



2 \a + h + + ah — 1 ah + a~h~ — 1 
Summarizing the computations thus far, one finds 

■ln(0+J = g+ -— In 



(3.52) 



<iz 2 <iz \ 9 J 2 \ z — b z — 6+ / z — 6+ z — 6 

' ' a~h~ ah z — a~hZah aha + hf — ah z a + h^ 



ah + a h — 1 a+/i + + ah — 1 

1 / a + ht — ah z ah z — a~h~ 



2 \a + h + + ah — 1 ah + a h — 1 / 
g+(S + -I)<P (3.53) 
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for some function $ = <£>(z, n). Using 

oc 

g(z,n) = G(z,n,n) = -^fji^z^' 1 , {z,n)eCxZ (3.54) 
j=o 

for l/\z\ sufficiently small (see, e.g., [221 Appendix C]) and l(3.32|l . one concludes 

oo oo oo 

- x>+^ j : = - e h z ~ j ^ + ( s+ j ) e ( 3 - 55 ) 

J=l J=0 j=0 

for some sequences {d+j ;(n)}„gZ) j G No- 
Noting that 



(3.56) 



one finds the analogous result regarding p_j and d_ j, j 6N. □ 
Remark 3.4. Closely related to (|3.53|) is the identity 



- In (l - - j = -2. 9 + (S+ - 7) - 2$ + - ln( 5 ) j , (3.57) 

which follows from l|3.44|l and (|3.53|l . see also equation (3.f9) in Equation 
(|3.57|l is t/ie analog of a well-known identity for the diagonal Green's function of 
Schrddinger operators, see, for instance, p. 369], |19| . |21l p. 99 and 122], and 

EH 

For later use we add the notion of degree of the quantities introduced. 
Remark 3.5. Introduce 

deg(a (r) ) = deg(6 (r) ) = 1, r G Z. (3.58) 

This results in 

deg(f t )=£, deg(&)=* + l, leN, (3.59) 

using induction in the linear recursion relations l|2.4f) - (|2.6|l . Similarly, the recursion 
relations i|3.27[) and 13.28(1 yield inductively that 

deg(0±, i±1 )=j, j£N . (3.60) 

4. Local conservation laws 

Throughout this section (and with the only major exception of Theorem 14. 4f> we 
make the following assumption: 

Hypothesis 4.1. Suppose that a,i:ZxM^C satisfy 

sup (\a(n,t p )\ + \b(n,t p )\) < oo, 
(n,t p )eNxR / 41 \ 

o(n, •)> ') S C^R), n G Z, a(n,i p ) ^ /or a/Z (n,t p ) eZxi 

To simplify notations we denote the bounded difference operator defined on 
£ 2 (Z), generated by the finite difference expression P2 P +2 in l|2.13|l . by the same 
symbol Pi P +i- 

We start with the following existence result. 
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Theorem 4.2. Let (z, t p ) € (C \ spec(iJ)) x K and suppose a, b satisfy Hypothesis 
\4-l\ and Tl p (a,6) = for some p £ No. Then there exist Weyl-Titchmarsh-type 
solutions ip± = ip±(z,n,t p ) such that 

V>±0, -,t p ) G f([n ,±oo) HZ), n G Z, ip±{z, n, • ) G C 1 (R), (4.2) 
and V± simultaneously satisfy the following two equations 

L(t p )ip±(z, ■ ,t p ) = zip±(z, - ,t p ), (4.3) 

i/)±,t p (z, - ,t p ) = P2 P +2(tp)i>±(z, - ,t p ), (4.4) 

= 2a(t p )Fp(z, ■ ,t p )ip^(z, ■ ,t p ) +G p+1 (z, ■ ,t p )ip±(z, ■ ,t p ). (4.5) 
Moreover, the Wronskian 

W(i()—(z,n,tp),il)+(z,n,tp)) is independent of (n,t p ) G Z x K. (4.6) 

Proof. Applying (H(t) - ziy 1 to <5o (cf. l|3.10j) and (^.lljl ) yields the existence of 
Weyl-Titchmarsh-type solutions ^>± of Lip = zip satisfying (|4.2(l . Next, using the 
Lax commutator equation (|2.26|) one computes 

z^±,t p = (m±) tp = L t J>± + LV± <t = [P 2p +2, i]*± + L^ ±t t 



zP2p+2^± ~ LP 2p+ 2^± + L^± ttp 



(4.7) 



and hence 

(L-zI){V±, tp -P2p + 2*±) = 0. (4.8) 

Thus, ^>± satisfy 

V±,t p - P 2p +2*± = C±*± + D ± ^ T . (4.9) 
Introducing *5f± = c±ip±, and choosing c± such that c±j — C±c±, one obtains 

i>±,t p - P2 P +2i>± = D±i/j t . (4.10) 

Since ip± G £ 2 ([n 0l ±oo) HZ), n G Z, and a, 6 satisfy Hypothesis 14.11 (|2.24H shows 
that P 2p +2V'± = {2aF p S + ip± + G p+1 ip±) G £ 2 ([n ,±oo) n Z). (Incidentally, this 
argument yields of course H4.5J) .) Moreover, since ip±(z,n 7 t p ) — d±(t p )(H(t p ) — 
zL)~ 1 5o)(n) for n G [±1, oo) O Z and some d± G C 1 (R), the calculation 

i>±,t p = d±,t p (H - zl)- x 8 - d ± (H - zL)- x U tv {H - z/)" 1 ^ (4.11) 

also yields tp±,t p G ^ 2 ([n ,±oo) n Z). But then D± = in (|4.10|> since ip T £ 
£ 2 {[n , ±oo) n Z). This proves JOJ. 

Since ip±(z, ■ ) satisfy Lip(z, ■ ) = zip(z, ■ ), the Wronskian (|4.(i|) is n- independent. 
To show also its i p -independence, a computation reveals 

■±-W(tl>-,1>+) = (°^ + 2{z-b + )F+ + G+ +1 + G p+ i)W-,V>+)=0 (4.12) 

using E3J), and (|2~35|) . □ 

For the remainder of this section, ip± will always refer to the Weyl-Titchmarsh 
solutions introduced in Theorem 14. 21 

The next result recalls the existence of a propagator W p associated with P 2p +2 ■ 

Theorem 4.3. Assume Hypothesis \4-l\ and suppose a,b satisfy Tl p (a,6) = for 
some p G No- Then there is a propagator W p (s, t) G £>(^ 2 (Z)), (s, t) G R 2 , satisfying 

(i) W p (t,t) = L, t€R, (4.13) 

(it) W p (r,s)Wp(s,t) =W p (r,t), (r,M)eR 3 , (4.14) 
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(Hi) W p (s,t) is jointly strongly continuous in (s,t) G R 2 , (4-15) 
such that for fixed to £ 1, /o 6 ^ 2 (^), 

f(t)=W p (t,t )f , 468, (4.16) 

satisfies 

j t f{t) = P 2p+2 {t)f{t), /(«„)=/„. (4.17) 

Moreover, H(t) is similar to H(s) for all (s,t) G R 2 , 

H(s) = Wp(s,t)H(t)W p (s,t)- 1 , (si)el 2 . (4.18) 

This extends to appropriate functions of H(t) and so, in particular, to its resolvent 
(H(t) — zl)^ 1 , z£C \ a(H(t)), and hence also yields 

a(H(s))=a(H(t)), (s,t)eR 2 . (4.19) 

Proof. (14.1 3(1 - (|4.17|l are standard results which follow, for instance, from Theorem 
X.69 of |28j (under even weaker hypotheses on a, b). In particular, the propagator 
W p admits the norm convergent Dyson series 

W p (s,t) = I + I [ dt 2 --- ( dt k P 2 p +2 (t 1 )P 2 p +2 (t 2 )---P 2 p +2 (t k ), 
keN Js Js Js 

(s,t)eR 2 . (4.20) 

Fixing s G R and introducing the operator-valued function 

K(t) = W p {s,t)H{t)W p (s, t)-\ teR, (4.21) 

one computes 

K'(t)f = W p (s,t)(H'(t) - [P 2p+2 (t),H(t)])W p (s,t)- x f = 0, t e R, / e £ 2 (Z), 

(4.22) 

using the Lax commutator equation l(2.26|l . Thus, K is independent of t G R and 
K(s) = H(s) then proves lj3"T%)l , □ 

In the special case where H(t), t G R, is self-adjoint, the operator P 2p+2 (t) is 
skew-adjoint, P 2p + 2 (t)* = —P 2p+2 (t), t € R, and hence W p (s,t) is unitary for all 

(M)el 2 . 

Next we briefly recall the Toda initial value problem in a setting convenient for 
our purpose. 

Theorem 4.4. Let t 0)P e R and suppose a(°\M°) e £°°(Z). Then the pth Toda 
lattice initial value problem 

Tl p (a,6) = 0, (a,b)\ tp=top = (a^,b^) (4.23) 

for some p G No has a unique, local, and smooth solution in time, that is, there 
exists a To > such that 

a(- ),&(•) G C 00 ((to,p-To,^ P + T ),£ 00 (Z)). (4.24) 

Proof. This follows from standard results in ^ Sect. 4.1] and has been exploited in 
the self-adjoint case in fjTI] Sect. 12.2]. More precisely, local existence and unique- 
ness as well as smoothness of the solution of the initial value problem (|4.23(l (cf. 
also follows from fIJ Theorem 4.1.5] since f p+1 and g p+ i depend only poly- 

nomially on a, b and certain of their shifts, and the fact that the Toda flows are 
autonomous. □ 
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Remark 4.5. (i) Theorem 4-1-5 in also shows that Theorem \4-4\ extends to the 
case where a 1 -- ^,^ ^ G £°°(Z) is replaced, for instance, by 

{« {0) (n?-\} n& ,{bW(n)} nez ef(Z). (4.25) 

To apply the Banach space setting of ^ Theorem 4.1.5] in this situation, one in- 
troduces a new sequence 

{c {0) W} n ^ = {" i0) W 2 -\} n ^ (4-26) 

and then demands 

{6 (0) (»)} ngz> {c (0) (n)} neZ e^(Z), (4.27) 

upon substituting -[o^(n) 2 } ngZ by {c^(n) + i} ngZ * n ^ e equations of the Toda 
hierarchy. The latter is possible by Remark \2-4\ since only a 2 but not a itself enters 
the Toda hierarchy of evolution equations. This observation will be used in the proof 
of Theorem \5. 7| 

(ii) In the special case where H(t), tGR, is self-adjoint, one obtains 

sup (\a(n,t p )\ + \b(n,t p )\) <2\\H(t p )\\ = 2\\H(t 0iP )\\, (4.28) 

(n,t p )£Nx(t a , p -T ,t 0tP +T ) 

using (the local version of) the isospectral property (|4.19|l of Toda flows. A further 
application of ^ Proposition 4.1.22] then yields a unique, global, and smooth solu- 
tion of the pth Toda lattice initial value problem (|4.23|) . Moreover, the same argu- 
ment shows that if a, b satisfy Hypothesis \4.1\ and the pth Toda equation T\ p (a, b) — 
0, then a, b are actually smooth with respect to t p G M., that is, 

a{n, ■ ), b(n, ■ ) G C*°°(R), n E Z. (4.29) 

Theorem 4.6. Assume Hypothesis \4-l\ and suppose a,b satisfy Tl p (a,6) = for 
some p e No. Then the following infinite sequence of local conservation laws holds, 

d tp ln(a) - (S+ - I)f p+1 = 0, (4.30) 

d tp P±, 3 + (S+ - I)J ±tPij =0, j G N, (4.31) 

where 

J±,P,j = -2 (a 2 M±,p+j-t) , (4-32) 

and 4>±,j, j G N, are given by <|3.27[l and (|3.28l) . and p±,j, j G N 7 are given by 
and (ESJ). 

Proof. Using (|4.5|l and (|3.21l) one computes 



ft„ In 



v4\ _ ^i,t, 



+ 

2a+F+0+ - 2aF p( j) + + G+ +1 - G p+1 



= (S+ - I) (2aF p + + Gp+i) . (4.33) 
Utilizing the expansion H3.26J) of 0+ as z — > oo one finds 

2aF p q>+ + G p+1 = 2a 2 U-lA ( £ t f f +d"~ i ) ~ zP+ ' + E Sr-t* 1 + fp+i 
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2a 2 ^2 fp-^+A + 9p + fp+i ) + ^2 ( 2fl2 E fp-fP+A + 9p-k J z k 



oo , p 



OO / p K 

— 1 + 2a 2 ( I] M+,P+j-i ) « 



(4.34) 



On the other hand, using l|3.21[l and the expansions l|3.26[) and <|3.32[1 . one concludes 



Combining these equations one infers 



p 

at 



3=1 v <?=1 

p 



- + Y^p+^ z ~ i = is+ - ^ 2fl2 E + 5p + /p+i 

^=1 

P , P X 

+ E ( 2fl2 E fp-i<t>+,t-k + g P -k ) z k 

k=l ^ £=fc+l 

00 p 

- + 2a 2 ]T]T , p+J _ £ z- 
j=i £=0 

= i s+ - !) ( 2a 2 ^2 fp-e<t>+A +9p + fp+i 



=1 
00 p 



+ (5+ - I) ( 2a 2 J] J] M+, P+] -tz- 3 ) (4.36) 
^ j=i e=o ' 

/ 00 p 

= (5+ - /)/ p+1 + (5+ - /) ( 2a 2 E M+a+i-t* 

^ j=i £=0 

and hence (|4.31(1 . 14.32|l in the case of p+.j- Here we used the first of the equations 
in Tip (a, 6) = as well as the fact that the left-hand side of (|4.36|) contains no 
positive powers of z. 

Similarly, one can start with 0_ and finds 

d tp ln f^f-) = (S+ - T)(2aF p ^ + G p+1 ) 



= ( s+ - I) ( 20-,o/ P + 2 ^2 ft<j>-,p-e +9 P + fp+i 
p-i , p-i-j 

+ (s+-i)j2h E + 20-,-i/p-j+i 

i=i ^ e=o 

+ 20-,o/p-j +5p-j^ j 
+ (S+ - I) ((2^_,_i/ p _i + 20_,o/o + .go)^ p + (20-,-i/o - l)^ p+1 ) 

OO , P s 

+ (5+ - /) J2 ( 2 E /^-.p+i-^ z ~ j - ( 4 - 37 ) 

i=i ^ £=0 ' 
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Thus, 



/ / 

+ »./' rr- ' = d tp ln(0 + ) = d tp In J 

j=l \ V 



+ 



= (5+ - I) (W-,o/ P + 2 ^ + 5 P + /p+i 

oo , p s 

+ (S+ - 7) £ ( 2 £ M-j+j-t ) z~ j , (4.38) 

implying 

dt p P-,j = (5+ - /) M-^ri-i) > 3 6 N (4.39) 

and hence (|4.31|) . H4.32|) in the case of □ 

Remark 4.7. W^e emphasize that the sequence (|4.39|l yields no new conservation 
laws, and so the latter is equivalent to that in l|4.36|l . 

The first local conservation law, 14.30fl . is of course nothing but the first equation 
in Tl p (a, b) — 0, namely a tp = a(f^ +1 — f p +i)- The second equation in Tl p (a, b) = 
is (|4.31|l for j = 1, namely 1 b^ = — (5+ — I)g p +i. The first few local conservation 
laws explicitly read as follows: 
(») p = 0: 

d toP+ , j = 2{S + -I)4> +d > (4.40) 

in particular, 

3 = 1: d to b+ = 2(5+ - I)a 2 , 

3 = 2: d to {{a + ) 2 + \{b+) 2 ^=2{S + -I)a 2 b+. 

(ii) p = 1: 

d tlP+tj = 2(5+ - I)(4>+. ]+ i + (b + ci)0 +J ), (4.41) 

in particular, 
3 = 1: d tl b + = 2(5+ - I){a 2 (b + + b)+a), 

J = 2 : d tl (V) 2 + i(6+) 2 ) = 2(5+ - /)(a 2 ((a+) 2 + b+(b + 6+)) + Cl 6+). 

Using Lemma \'3.tA one observes that one can replace p±j in l|4.31|l by fj by suitably 
adjusting the right-hand side. 

An obvious consequence of the local conservation laws is that, assuming sufficient 
decay of the sequences a 2 — j and b, one obtains 

J- p±j ( n ' t p) = 4rll h *p) = °- ( 4 - 42 ) 

P nGZ P n€Z 



^Using the second equation in Tl p (a, b) = one infers that 2a 2 5Z?=o flfi+yP+l—l = ~9p+l 
up to a constant, that is, an element in the kernel of S + — I. Using the notion of a degree (cf. 
Remark 13. 51 . one concludes that the constant equals zero. 
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5. HAMILTONIAN FORMALISM, VARIATIONAL DERIVATIVES 

We start this section by a short review of variational derivatives for discrete 
systems. Consider the functional 

T: £°°(Z) K -> C, 

for some k £ N and k £ No, where F: Z x C 2rK — > C is C 1 with respect to the 2rn 
complex- valued entries and where 

u^=S^u, SW = { {S+ J™ ifS " 0, ue£^(Zr. (5.2) 

I (S 1 ) s u if s < 0, 

For brevity we write 

F(n, u) = F(n, u(n), (n), (n), . . . , u (fe) (n), (n)) , (5.3) 

and it is assumed that {F(n,u)} n £Z £ £ 1 (Z). 

For any w £ ^ 1 (Z) K one computes for the differential dT 

{dT) u (v) = ^H u + eu)| e=0 

flF(n,u) (fc) gf(n,u) ( _ fc) 
+ + cM fc ) W+ du(~ k ) { ' 
( dF(n,u) ^ dFfau) (+1) dF{n,u) 



assuming 



{F(n, M )} neZ , { ^"'^ j ^(Z), j = l,...,fc. (5.5) 

Because of the result 1)5 .4JI . we thus introduce the gradient and the variational 
derivative of T by 



assuming 115.5(1 . 

To establish the connection with the Toda hierarchy we make the following as- 
sumption for the remainder of this section. 

Hypothesis 5.1. Suppose 

a, a' 1 , bee°°(Z). (5.7) 
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We introduce the difference expressions 

D= {p 2 D 1 =aS+ -a, D 2 = a-a-S-, 

D -i-( o (D , ( 



where 

n-1 



TYl— — OO 



my u )(n) = — E = E S «e^(Z). 

m— n+1 m— n 

Viewing D,D~ X and Di, , Z?2, -D^ 1 as operators on £ 1 (Z) 2 and ^(Z), respec- 
tiveiy, one conciudes that 

DD- 1 = 1*10)2, i^- 1 = / /MZ) , Z^I^ 1 = 7 £1(z) . (5.10) 

Next, let T be a functional of the type 

J 7 : r°(Z) 2 -»C, 

T{a, b) = E a, &, ■ ■ ■ , a (+fc) , a ( - fc) , 6 ( - fe) ) 

nGZ 
riSZ 

assuming 

(5.12) 

The gradient and symplectic gradient V S T of T are then defined by 

^ - ((v?!;) - (|) ^ 

and 

(V s ^) a , fc = D(VF) a , b = D , (5-14) 

respectively. In addition, we introduce the bilinear form 
ft: l\Z) 2 x ^(Z) 2 -» C, 

v) = lYl {( d ~ 1u ) (") • + • ((£> _1 )Mw) • (5.15) 

nez 

One then concludes that 

£l(Du, v) — E] u(n) ■ v(n) = E] ( u i( n ) v i( n ) + ii2(i)t>2(i)) 

nez nez (5.16) 
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where ( • , • denotes the "real" inner product in £ 2 (Z) 2 , that is, 

(■,-) P(z y.l 2 {1) 2 x£ 2 (Z) 2 ->C, 

(u,v)p(zy = ^2u(n) ■ v(n) = ^ (wi(n)ui(n) + « 2 (n)u 2 (ra)). ( 5 - 17 ) 

nGZ nSZ 

In addition, one obtains 

{dT)a,b{v) = {{VT) aM v) t 2 (z) 2 = n(U(VJ0« l6 , «) = n((V.J0«, 6) «)• (5.18) 
Given two functionals T\ , .7-2 we define their Poisson bracket by 
{T X ,T 2 } = dF x {V s F 2 ) = n(V a ^i,V,^ a ) 

- S tel.!) (*!-!)■ 

One then verifies that both the Jacobi identity 

{{J' 1 ,^ 2 },^ 3 } + {{J' 2 ,^ 3 },^ 1 } + {{J'3,^i},^ 2 } = 0, (5.20) 
as well as the Leibniz rule 

{^1,^3} = {^1,^ 2 }^3 + ^ 2 {^1,^3}, (5.21) 

hold. 

If is a smooth functional and (a,b) develops according to a Hamiltonian flow 
with Hamiltonian H, that is, 

fj =(V s W) a , h = #(VW) Q , b = #(g) , (5.22) 



<5b 



then 



n 



dt dt 

= k \%v>) ' w)A - h {% wJ ' 

= {^,W}. (5.23) 

Here, and in the remainder of this section and the next, time-dependent equations 
such as H5.23(l are viewed locally in time, that is, assumed to hold on some open 
i-interval ICR. 

If a functional Q is in involution with the Hamiltonian TL, that is, 

{5,H} = 0, (5.24) 

then it is conserved, that is, 

f = 0. ,5, 5 , 
Next, we turn to the specifics of the Toda hierarchy. 
Lemma 5.2. Assume Hvvothesis \2. 1\ Then, 

6 -p- = --gt-u ^H, (5.26) 
da a 
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5 -lj-=lf t -i, ten. (5.27) 

Proof. With our assumptions on (a,b) we only know that fg € We start by 

deriving (|5.27[1 . To that end we introduce the functional 

hN{a,b) = Y,h(n)XN{n), (5.28) 

nGZ 

where xn is the characteristic function of the set [-N, N] n Z. Then one finds 

at (a, 6) ) 6 (u) 

- E twf + xnW^Sv + ^ s(+1) wk + • • •) 

raGZ v 7 

= E^( n M n )> «e^(Z)- (5-29) 

nGZ 

On the other hand, recalling that (•, - )/3(z) denotes the usual scalar product in 

(dF i>N (a,b)) b (v) = -r^LN{b + ev)\ €=0 = V ( &n, L k vL t ~ 1 ~ k S n ) x»W 
de nGZ V fe=o '< 2 (z) 

= E f^'E^ 1 ^") Xjv(") + E f^^El^^ 1 ^]^) Xiv(n) 

nGZ V fc=0 J nGZ V fe=0 ^ 2 (Z) 

^ E E (*»> + E E ^"W^ 

fc=o nez fc=o nez 

= £^(5 n ,^- 1 5 n )^ (z) «(n) =^53/ < _i(n)«(n), « 6 ^(Z), (5.30) 

nGZ nGZ 

using 12.11|) . and the general result that for bounded operators A,B £ B(Tt) on 
a separable, complex Hilbert space 7i with AS and -BA trace class operators, 
their commutator is traceless, that is, tr([A, B]) — (cf. [IT]. [2*91 Corollary 3.8]). 
Combining the two expressions, one concludes that (|5.27|l holds. 

Next we turn to the proof of (|5.26|) . To this end one first observes that, as before, 

^ E (f + s m, #t> + s<«>^ + )»w 

nGZ v ' 

= Y^l(n)v(n), vtlHZ). (5.31) 
oa 

nGZ 

Furthermore, one computes 
(d!Pi iN (a,b)) a (v) = ^-?£,v(a + ew)| e=0 



^E (( a + e ' i; )'5' + + (a + eu) 5 + b) 1 8 n ) p ^xn (n) 

nGZ 



e=0 
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2 (L k (vS+ + v-S-)L e - 1 - k )S n ) X N 

t-x 

£-1 

+ E E (*«> [ ifc ( v5+ + u ~ 5 ~)> ^- 1-k ]«5»)xw(») 



fe=0 n6Z 
£-1 



^ E E + 

k=Q neZ 

+ ^tr([i fc (^+ + W -^-),L £ - 1 - fe ]) 

fc=0 

I («(n) (5-«„, ^-^„)^ (z) + «(n - 1) (S+<5„, i^^n), 2(z) ) 
^E ("W^+b^'^iip) l)(5„_i,i £_1 (5„ 
^E ( w ( n )( <S «+i' L ^ 1(5 ™)£2 (z) +u(ri)(5 n ,L^ _1 5 n+ i)^ 2(z) ^ 
2£^(<5„ +1 ,L £ - 1 ,5„) £2(z) ^) 



= -^-^j w (n)»(>i), WG^CZ). (5.32) 
t-^ a(n) 

Thus (l5~251) is established. □ 



For the remainder of this section we now introduce the following assumption. 
Hypothesis 5.3. In addition to Hypothesis \2. 1\ suppose that 

{a(n) 2 - l) nez ,{b(n)}nez e (5.33) 

We fixleN and want to show that a suitably renormalized functional Ti is 
well-defined. We define it by subtracting the limit of each ft(n) as n — > oo. Each 
fe is a polynomial in a, b and certain shifts thereof. Terms with a and shifts enter 
homogeneously and only in even powers. Assume that fe(n) — > A^ as |n| — > oo. 
For £ odd, we know that A^ = because each each term contains at least one b or 
certain shifts thereof. Then {fe{n) — Ag} ne z G ^(Z). We can see this as follows. 
Each fa is a finite sum of terms containing a and b with shifts. Terms with b 
are already summable. Only terms exclusively with a have a nonzero limit and 
hence are nonsummable. Their general form will be «i(n) • • • afe(n) with otj(n) = 
as \n\ — > oo, where (rny) denotes the shifts. Then 

ai(n) ■ ■ ■ afc(n) - Ai • • • A* = (ai(n) - Ai)a 2 (n) • • • otfc(n) 

+ Ai(a 2 (n) - A 2 )a 2 (n) • • • afc(n) 
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+ --- + A 1 .--A k _ 1 (a fc (n)-Afc). (5.34) 

Thus we see that 

{ai(n) • • • a k {n) - \ x ■ ■ • A fc }„ eZ G (5.35) 

and hence 

{/*(*)- Mnez G * X (Z). (5.36) 

Thus, the functional 

?t(a,b) = '%2(ft(n)-\t) (5.37) 

is well-defined. 

In the following it is convenient to introduce the abbreviation N_i = NU{ — 1, 0}. 

Theorem 5.4. Assume Hvvothesis \5.,'A and let p G N_i. In addition, define 

= ^ E (WH - V), (5-38) 
p i ^ 



ne- 



where 



Then, 



Ap+2 = , lim fp+2(n). (5.39) 

\n\ — >oo 



^-JTi^—Z^- (vfp,l = ?T2^ a = A + " < 5 ' 40 > 

Moreover, the homogeneous pth Toda equations then take on the form 

%{a,b)= (*>) -D[ {VI l] a ) =0, peN-j. (5.41) 



Proof. This follows directly from Lemma f5. 21 and ()2.29fl . □ 
Remark 5.5. in <|5.41[) we also introduced the trivial linear flow Tl_i given by 

a -' ( -«-te)- fl (S: 1 1 h)-te)-' 1 ,5 - 42) 



Next, we consider the general case, where 

't P - a(/+ +1 - ,/ 
+ 9 P +i ~ 9p+i 



Tl p (a,b)=( at l U+l) )=0, pen-,, (5.43) 



and 

v v 
fp = ^2 c p-kfk, 9p=y~] Cp-k9k - Cp+i, p G N . (5.44) 
fe=0 fc=i 

Introducing 

p 

/-!=/-!, I P = ^c p _ fc 4, peN , (5.45) 

one then obtains 

Tl p (a, b) = ( a A - D (Z?:) peN-,. (5.46) 



(yi P )b 
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Theorem 5.6. Assume Hvvothesis \5 .'A and suppose that a,b satisfy Tl p (a,6) = 
for some p € N_i. Then, 

" Ir = 0, re N_i. (5.47) 



dt p 



Proof. From (|4.31l) and 1|3.39[) (cf. Remark l4~7jl one obtains 

{S+-I)P r+2 , reN-i, (5.48) 



dfr+2 



dtp 

for some (3 r +2, r G N-i, which are polynomials in a and 6 and certain shifts thereof. 
Using definition (|5.38[1 of I r , the result (|5.47|) follows in the homogeneous case and 
then by linearity in the general case. □ 

Theorem 5.7. Assume Hvvothesis \5.3\ and letp,r G N_i- Then, 

{I p ,I r } = 0, (5.49) 
that is, I p and I r are in involution for all p,r G N-i. 

Proof. By R,emark l4.5l (i). there exists To > such that the initial value problem 

T\ p (a,b) = 0, (a, 6)| t=Q = (a' '^), (5.50) 

where , satisfy Hvpothesis l5.3l has a unique and continuous solution a(t), b(t) 
satisfying Hypothesis 15 . 31 for each t G [0, T). For this solution we know that 

j t I p (t) = {I r (t),I p (t)} = 0. (5.51) 

Next, let t I 0. Then 

= {I r (t),Ip(t)} ^ {/ P (0),/ p (0)} = {Ir,I P }\ [atb)=(amtbmy (5-52) 

Since a^°> , are arbitrary coefficients satisfying Hypothesis 15.31 one concludes 
(EH. □ 

There is also a second Hamiltonian structure for the Toda hierarchy. Rewriting 
the linear recursion $n $ -$T ty in the form 

ft+i - ft+i = -\iSt - 9j) + (b + ft - bh), 

g (+1 -g w = -2(a 2 f+ ~(a-) 2 f^) + b( gi ~gi), £ G No, (5.53) 

one finds 

2"Un+l - /n+l)\ _ jj f^a9n\ ^\ 
9n+l-9n+l J V fn J ' 



where we abbreviated 



§a(S+-S-)a a(5+-J)6 
6(J-S-)a 2(S'+-S'-)a 2 



One can introduce a second Poisson bracket {{ • , • }}, defined by 



(5.56) 
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This second Poisson bracket is also skew-symmetric and satisfies the Jacobi iden- 
tity and the Leibniz rule. As in Theorem 15. 71 one verifies that all I p , p € N_i, are 
in involution also with respect to this second Poisson bracket (|5.56|) . that is, 

{{I r ,I p }} = 0, p.reN-r. (5.57) 

Combining Ij5.38jl . (|5.45(l . (|5.46|) . and H5.53J) . 15.55|l permits one to write the pth 
Toda equation in another Hamiltonian form 



6. A SKETCH OF THE ALMOST PERIODIC CASE 

Following Johnson and Moser |2U in the context of one-dimensional almost pe- 
riodic Schrodinger operators and Carmona and Kotani ^ in the corresponding dis- 
crete case (see also the treatment in Carmona and Lacroix Sects. VII. 1, VII. 2]), 
we now very briefly sketch an extension of the Hamiltonian formalism to the case 
of almost periodic coefficients a and b. 

To set the stage, if c denotes an almost periodic sequence on Z, we recall that 
its ergodic mean (c) is given by 

JV 



c) = lim c„. (6.1) 

1 vtoo 2N+ 1 ™ v ; 



n=-N 



We assume that u has the frequency module M.(u) and given a density F as in 
(15.31) . equation (|5.1I) then becomes 



N 

F(u)=Yim F(n,u(n),u( +1 \n),u ( >- 1 '>(n),...,uW(n),u ( -- k '>(n)), ^ 

and assuming that the frequency module A4(v) of v satisfies M.(v) C A4(u), the 
analog of (|5.4|) then reads 

{dT) u (v) = -^^"(u + ™)| £=0 

N 



1} dF(n,u) 

-N 



lim Y ( d J^ + S ^ d 4^ + S^ 



+ 



= ((Vn,.»> = (^»). (6.3) 

Assuming a and b are almost periodic with frequency module M., we again 
consider functionals as in (|5.11|) . 

F{a, b) = (F( ■ , a, b, a( +1 \ . . . , a^ +k \ b^ +k \ a^ k \ &<-*>)) 

= (F(-,a,b)), (6.4) 
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and in analogy to (|5.19(l , the Poisson brackets of two functionals T\ , Ti are then 
given by 

<*-MK); fl (*9>- 

Again one verifies that both the Jacobi identity as well as the Leibniz rule hold in 
this case. Moreover, as in (|5.22|) and l|5.23[l . if T is a smooth functional and (a, 6) 
develops according to a Hamiltonian flow with Hamiltonian TL, that is, 

ty= (V s W) a , fc = D{VH) a , b = D f ||j , (6.6) 

then 

f = |<F(n,M)> = {W (6.7) 

Next, assuming in addition that 1/a 6 £°°(N) and that (1/a) stays finite for 
t varying in some open time interval 1 C 1, we now introduce the fundamental 
function w by 

w(z) = (4> + (z, •)} (6.8) 

for \z\ sufficiently large. Since 

w'{z) = -(g(z, ■)), zeC\spec(tf), (6.9) 

w in (|6.8(l extends analytically to z G C \ spec(iJ). We note that w is the ap- 
propriate discrete analog of the function (m + {z)) = (— 1/G(z, • , • )) introduced 
by Johnson and Moser [22] in the case of almost periodic Schrodinger operators. 
(Here m + denotes the right half-line Wcyl-Titchmarsh coefficient and G(z, • , • ) the 
diagonal Green's function of the underlying almost periodic Schrodinger operator 
H = -d 2 /dx 2 + V in L 2 (R).) 

One observes the asymptotic expansion 



(6.10) 



w(z) = (ln(cj> + (z,-)) = 0n(a/z))+5> +J >ar- 
= (ln( a /^))+^ J - 1 </,> 

z— >oo * * 

jen 

and computes 

-^^-(n) = — V — - 2G(z,n,n+ 1), z G C \ spec(iT), n G Z, (6.11) 
oa a(n) 

^jp(n) = -p(z,n), zeC\spec(ff), neZ. (6.12) 
Introducing 

1 p 
I-i=I-i= (A), J p = — -(/ p+2 >, I p = £cp_ fe / fc , peN , (6.13) 

p + ^ fe=0 
the Toda equations again take on the form 

tpj \(v/ P ) b/ 



Tl p (a,b) = )-D[) T7r %" )=0, pGN-i, (6.14) 
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Finally, one can show that w(z\) and w{z2) are in involution for arbitrary zi, zi £ 
C \ spcc(iJ), and hence obtains 

{w(zi),w(z 2 )} = 0, zi,z 2 e C\spec(fl"), (6.15) 

= 0, {I p ,I r } = 0, zeC\spec(ff),p,reN_i. (6.16) 

Naturally, these considerations apply to the special periodic case in which (c) 
for a periodic sequence c on Z is to be interpreted as the periodic mean value. 
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